Dwarf Galaxy Formation Was Suppressed By Cosmic Reionization 
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A large number of faint galaxies, born less than a billion years after the big 
bang, have recently been discovered The fluctuations in the distribution of 
these galaxies contributed to a scatter in the ionization fraction of cosmic hydro- 
gen on scales of tens of Mpc, as observed along the lines of sight to the earliest 
known quasars Theoretical simulations predict that the formation of dwarf 
galaxies should have been suppressed after cosmic hydrogen was reionized i'^, 
leading to a drop in the cosmic star formation rate Here we present evidence 
for this suppression. We show that the post-reionization galaxies which pro- 
duced most of the ionizing radiation at a redshift z ~ 5.5, must have had a mass 
in excess of ~ lO^^^^^'^Mo or else the aforementioned scatter would have been 
smaller than observed. This limiting mass is two orders of magnitude larger than 
the galaxy mass that is thought to have dominated the reionization of cosmic 
hydrogen (~ lO^Mo). We predict that future surveys with space-based infrared 
telescopes will detect a population of smaller galaxies that reionized the Universe 
at an earlier time, prior to the epoch of dwarf galaxy suppression. 

The lines of sight towards the earliest known quasars at redshift s 2; > 6 show a rapid 
transition in the ionization state of the intergalactic medium (IGM), potentially marking 
the end of the reionization epoch ^^-^^ , It is thought that the dominant source population 
to have triggered reionization included dwarf galaxies with virial temperatures above the 
hydrogen cooling threshold (Tvir ~ lO^K). The re-ionization of cosmic hydrogen resulted in 
heating of the IGM to ~ lO^K which drastically increased the minimum virial temperature 
of new galaxies ^^"^'^ to a value as high as Tvir ~ lO^K. Reionization is therefore expected 
to be followed by a depletion of low-mass galaxies and hence a decrease in the global star 
formation rate The detection of this anticipated suppression at 2; ~ 6 would provide 
important evidence for a late reionization epoch. 

In addition to luminous quasars, starburst galaxies with star formation rates in excess of 
~ O.IM0 yr~^ and dark matter halos ^ of ~ 1O^~^°M0 have recently been discovered 1-6,19-22 
at 2; ~ 5-6. These sources contributed to the ionizing background radiation at the end of 
the reionization epoch. Luminous Lya emitters are routinely identified through continuum 
dropout and narrow band imaging techniques ^,4,20 However, in order to study fainter 
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sources which were potentially responsible for reionization, spectroscopic searches have been 
undertaken near the critical curves of lensing galaxy clusters 21-23^ where gravitational mag- 
nification enhances the fiux sensitivity. A comparison between the findings of these deep 
surveys and the statistics of bright Lja emitters provides a preliminary hint for the deple- 
tion of low mass galaxies, as expected from galaxy formation in a photo-ionized IGM i2-i3_ 
However, this inference relies on a comparison between heterogeneous data sets based on 
different methods of discovery. 

The upper right panel of Figure 1 shows data on the luminosity function of galaxies at 
z ~5.5-6 from the Hubble Space Telescope Ultra-Deep Field (UDF) We first use this data 
to constrain the physical characteristics of high redsliift galaxies (see § 1 of the Supplemen- 
tary Online Materials) . We model starburst galaxies based on a simple prescription for their 
continuum luminosity combined with the Sheth-Tormen mass-function of galaxy halos. 
The two free parameters in this model are the starburst lifetime (tit) and the mass fraction of 
virialized baryons which are converted into stars in each galaxy (/star)- The lower right panel 
shows contours of the joint a-postcriori probability distribution [d"^ P / d(log fstm )d{log tit)] (see 
§ 2 of the Supplementary Online Materials for a description). The best-fit model luminosity 
function is plotted over the data, with the favored values of tit and /star listed. The other con- 
tours (grey) show the virial temperature corresponding to the lowest luminosity bin [labeled 
by logio(Tvir)]. The best-fit luminosity function imphes that the faintest observed galaxies 
have a virial temperature ~ lO^K, and that starburst galaxies at 2; ~ 6 have duty-cycles 
of tit/tu ~ 10% (where tu is the Hubble time), combined with star formation efficiencies of 
/star ~ 10 — 15%. Interestingly, a duty-cycle of tn/tu ~ 0.1 corresponds to the dynamical 
time at the virial overdensity of 200, i.e. the time it takes the virialized gas to settle to 
the center of a collapsing galaxy. For the purpose of understanding the process of galaxy 
formation at high redshifts and the way it shapes reionization, it is crucial to know whether 
star formation was suppressed in dwarf galaxies near the hydrogen cooling threshold. With 
a virial temperature of Tyir ~ lO^K, these galaxies had a total (halo) mass of ~ IO^Mq at 
2; ~ 6. Unfortunately, present detection limits are insufficient to probe the existence of such 
galaxies. 

We overcome the limitation of direct detection of faint galaxies by making use of the 
fact that the absorption spectra of high redshift quasars reveal large fluctuations in the IGM 
ionization state on very large scales '^-^'^'^ (tens to hundreds of Mpc). While observations of 
the galaxy luminosity function are flux limited, the fluctuations in the ionizing background 
are sensitive to contributions from galaxies of all luminosities. A measurement of the fluc- 
tuations in the ionizing background therefore provides a powerful probe of the properties 
of the galaxies generating the ionizing background radiation at the end of the reionization 
era, even if those galaxies themselves lie below current detection limits. The background 
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fluctuations can be inferred from observed fluctuations in the effective Lya optical depth, 

Tcfj. Fan ct al. ^ have measured average values for the effective optical depths of = 2.5, 
2.6, 3.2 and 4.0, within redshift bins of width Az = 0.15 at 2; = 5.25, 5.45 5.65 and 5.85, re- 
spectively. At these redshifts the scatter in t^s among the different lines-of-sight is a,- = 0.5, 
0.6, 0.8 and 0.8, respectively. The dominant uncertainty in individual estimates of Tcs is 
the level of unabsorbed quasar continuum, which introduces noise ^ in Teg of ±0.05. This 
observational uncertainty must be subtracted (in quadrature) from a-,- to obtain an esti- 
mate of the intrinsic scatter in Tes- Following this operation we find fractional fluctuations 
Sr = (Tr/ns of (0.20 ± 0.03), (0.23 ± 0.03), (0.25 ± 0.03), and (0.20 ± 0.03) at the above 
sequence of redshifts. The corresponding comoving mean-free-paths of ionizing photons are 
i^mfp = (39ll^)Mpc, (4ll2i)Mpc, (27l^^)Mpc and (19lig)Mpc respectively ^. Here we have 
corrected the frequency averaged ionization cross-section for a spectrum corresponding to a 
star-forming galaxy At the highest redshift under consideration, some lines-of-sight con- 
tain blank absorption troughs making their estimated value of Tcs a lower limit. At 2; ~ 5.85 
there is no measured lower limit for -Rmfp, and so we conservatively set it to 3Mpc (below 
the typical value at 2; ~ 6). 

We have calculated the predicted scatter in from a combination of cosmic variance 
among different regions and Poisson noise in the number of ionizing sources, based on an 
extension of a detailed model pubhshed elsewhere A description of this model is provided 
in § 3 of the Supplementary Online Materials, and only a summary of its main features 
is given here. Fluctuations in both the ionizing background intensity, J, and the density 
field around their mean values, contribute to fluctuations in Teff. Within our model, the 
ionizing background is smoothed on the scale -Rmfp, since any given point in the IGM is 
exposed to all sources within a volume V^fp = (47r/3)-R^jjfp around it. The model combines 
four major sources of fluctuations: (i) Poisson noise in the number of sources contained 
within a volume V"mfp; (H) Fluctuations in the value of J within vohimcs of V^mfp resulting 
from delayed (or enhanced) structure formation in underdense (overdense) regions In 
evaluating this quantity we include two contributions to the derivative dJ/dz arising from: 
(a) the evolution of the star formation rate, which at first we take to be proportional to 
the derivative of the mass fraction of baryons that were assembled into galaxies; and (b) the 
evolution of -Rmfp! (Hi) Fluctuations in Lya transmission due to fiuctuations in the density 
contrast smoothed on the scale -Rmfp; (iv) Fluctuations in the transmission due to small-scale 
density fluctuations along the line-of-sight through a region with constant J. These latter 
fluctuations are computed based on the results of numerical simulations 

We have compared the predictions of the above model with observational data as a 
function of the lifetime tit and the minimum virial temperature Tmm of the lowest mass 
galaxies that make a major contribution to the ionizing background. Contours of the re- 
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sulting distribution d^P/d(logtit)d(logTniin) (see § 4 of the Supplementary Online Materials 
for details) evaluated at z ~ 5.45 (blue contours) and z ~ 5.65 (red contours) are plotted 
in the lower left panel of Figure 1. The corresponding marginalized differential probability 
distributions rfP/c? ( log Tmin) are plotted in the upper left panel (solid blue and red lines). 
In addition we plot (iP/(i(logTmin) at 2; ~ 5.25 (dashed blue line) and z ~ 5.85 (dashed 
red line). Intriguingly, the results ai z < 5.65 favor Tmm ^ lO^K, while results at 2; ~ 5.85 
favor Tmin ~ 10^~^K [though note that (7^(2; = 5.85) is a lower hmit]. This trend is expected 
if reionization completed at 2; ~ 6-7. Indeed, the mini-halos (Tyir < lO^K) which provide 
the collapsed gas reservoir for star formation in galaxies above the hydrogen cooling thresh- 
old are thought to be photo-evaporated on a timescale of ~ OMu, which corresponds to 
Az ~ 1. Thus, if reionization completed at 2; ~ 6.5 — 7, we would only expect to observe the 
strong effect of dwarf galaxy suppression at 2; < 5.5 — 6. We also plot the combined differ- 
ential and cumulative distributions using data from the three redshift bins below z ~ 5.75 
(black lines), and find values of T^in ~ lO^-SiO-Sx (68%), with Tmin < lO^K ruled out at a 
(statistical) confidence level of > 99%. The above constraints have conservatively assumed a 
prior probability distribution for T^^^ that extends down to Tmin = lO^K. On the other hand, 
the ratio between the probabilities of having lO^K < Tmin < lO^K and lO^K < T^i^ < lO^K 
is independent of any assumed lower cutoff. We find this ratio to be ~ 0.05. A limiting 
virial temperature of Tmin > lO^K is an order of magnitude above the threshold for star 
formation triggered by atomic hydrogen cooling. However, Tmin ^ lO^K is consistent with 
galaxy formation in an IGM that was heated by UV photo-ionization at a higher redshift, 
thus preventing starbursts in low-mass galaxies from dominating the global star formation 
rate hy z 5.5. 

Our analysis so far has assumed for simplicity that the star formation rate is proportional 
to the derivative of the collapsed fraction in halos with Tyir > Tmin- However, starburst 
episodes arc believed to be triggered by galaxy mergers. We have therefore repeated the 
above calculation, replacing the derivative of the collapsed fraction with a star formation 
model that combines the contribution from mergers of halos with an initial Tyir < Tmin 
(that go above Tmin after the merger) and the contribution from mergers of galaxies with an 
initial Tyij. > Tmin (see § 5 of the Supplementary Online Materials for details). The results 
are shown in Figure 2 of the Supplementary Online Materials. This alternative estimate of 
the star formation rate also leads to a preferred value of Tmin ~ lO^'^^^'^K, with Tmin ^ lO^K 
at the > 99% confidence level. 

In the local Universe, star formation within galaxies with stellar masses below ~ 10^°Mq 
is suppressed as a result of supernova-driven winds which expel gas from their shallow po- 
tential wells high redshifts the lifetime of massive stars could be comparable to 
the halo dynamical time and this feedback may not operate as effectively. Nevertheless, it 



- 5 - 



is important explore the effect of such a feedback on fluctuations in the ionizing background 
and hence on our conclusions regarding the suppression of dwarf galaxy formation in the 
reionizcd IGM. We have repeated our analysis of the galaxy luminosity function, with an 
allowance for a suppression of star formation in galaxies below a critical mass (which is 
treated as a free parameter). Details of the required modifications and the results are de- 
scribed in § 6 and Figure 3 of the Supplementary Online Materials. Allowing for feedback in 
high redshift galaxies leads to estimates of /star ~ 0.5, with a duty-cycle that is larger than 
tit/tu ~ 0.5. The best fit value for the critical mass scale below which supernovae feedback 
operates is ~ 2 x lO^^M©, or TVir ~ lO^K. 

Estimates of the scatter in the ionizing background due to Poisson noise in the number of 
galaxies and due to fluctuations in the star formation rate as a result of delayed or enhanced 
structure formation must also be modified to account for possible supernova feedback. We 
find that if supernova feedback operates at high redshift in the same way as observed at 
low redshift, then the most likely value is Tmin ~ lO^'^-o-^K in agreement with our previous 
calculation, while Tmin ^ lO^'K is ruled out at the > 95% level. Note that the critical mass- 
scales derived for supernovae feedback and Mmin are similar, although they were derived 
from difi^erent data sets and astrophysical considerations. However, our results imply that a 
sharp cutoff due to reionization is required to explain the observed scatter in the effective 
optical depth of the IGM in addition to any suppression resulting from supernovae feedback. 

The observed star formation rate in the post-reionization era is on the borderline of 
being sufficient to reionize the Universe The value of Tmin ^ lO^K at z ^ 5.5 happens 
to describe the faintest galaxies in the UDF. We therefore predict that star formation is 
suppressed in galaxies just below the current limit of galaxy surveys at 2; ~ 6. Future 
searches at higher redshifts with JWST or large-aperture ground-based infrared telescopes 
should find a population of galaxies with T^ir ~ lO^K and detect an enhanced star formation 
rate prior to completion of reionization 
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Fig. 1. — Constraints on the star-formation efficiency (/star), duty-cycle (tit/tn) and mini- 
mum virial temperature (Tmin) of galaxies at z ~ 5.5 — 6. Panel a: The luminosity function 
inferred from the UDF together with the best-fit solution for a luminosity function model 
based on the Sheth-Tormen mass function (with luminosity proportional to halo mass). 
■M.1350 is the galaxy's absolute AB magnitude at a rest-frame wavelength of 1350A. The 
upper axis shows the halo mass corresponding to Mis^o for the best fit. Panel b: Con- 
tours (at 64, 26, and 14% of the maximum likehhood) of the joint a-posteriori probability 
distribution [(PP/d{log fstax)d{^ogtit)]. For comparison, the grey contours show the virial 
temperature of halos [labeled using logiQ(Tvir)] corresponding the lowest luminosity bin in 
the UDF. Panel c: The joint probability distributions for log(tit/tH) and log(Tinin), which 
include the constraint on tit/tn from the galaxy luminosity function. The contours represent 
values at 64, 26 and 14% of the maximum likelihood, and are shown for observations at 
z ~ 5.45 (blue contours) and z ~ 5.65 (red contours). This distribution includes constraints 
on tit from the galaxy luminosity function. However assuming the most likely value of -Rmfp, 
we find that Tmin ^ lO^K at z < 5.5 in all cases where the duty-cycle of the UV-bright phase 
of galaxies is greater than 1% (corresponding to the lifetime of massive stars). Panel d: 
Differential probability distributions for log (Tmin) at z ~ 5.25 (dashed blue lines), z ~ 5.45 
(blue lines), z ~ 5.65 (red lines) and z ~ 5.85 (dashed red lines). We also show combined 
constraints using the lowest three redshift bins (black lines). For comparison, the upper 
axis shows the corresponding masses sX z = 5.5. In this case the solid and dotted lines 
correspond to differential and cumulative distributions. Throughout this Letter, we adopt 
the latest values for the cosmological parameters as inferred from the Wilkinson Microwave 
Anisotropy Probe data 
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SUPPLEMENTARY ONLINE MATERIAL 

The above Letter to Nature interprets the observed luminosity function of z ~ 6 galaxies 
and the large-scale fluctuations in the Lya optical depth towards z ~ 6 quasars, based on a 
detailed theoretical model. Below we provide full details on the basic version of this model 
as well as its variants. The following sections expand the brief descriptions that can be found 
in the original Letter to Nature. 



1. Modeling the galcLxy luminosity function 

The properties of the observed starburst galaxies are interpreted through a comparison 
with a model luminosity function. We employ two estimates of the relationship between 
the luminosity of a starburst and its host halo mass. Our fiducial model assumes that 
the starburst luminosity L is proportional to halo mass M, i.e. L oc M. In this case, the 
luminosity function ^ (with units of number of galaxies per magnitude per comoving volume) 
may be written as 

^ tit dn dM 
~tndMdMis5o' 

where Ali35o(-^) is the starburst luminosity expressed in AB magnitudes at a rest- frame 
wavelength of 1350 A. Following Loeb et al. (2005), we express the relation between lumi- 
nosity and halo mass through two free parameters: (i) the mass fraction of virialized baryons 
converted into stars (the so-called "star formation efficiency"), /star; and (ii) the fraction of 
the Hubble time during which the starburst emits most of its ionizing radiation, tit/tn (the 
so-called "duty-cycle"). 



2. A-posteriori constraints from the luminosity function 

Using the models described above, the joint a-posteriori probability distribution ^^f^^^^^ 

for Usr and t^t was calculated as oc ^/star,t,t^^^^> where Lf^^^^^tn is the likehhood 

for the combination of /star and tit given the observed luminosity function. Here we have 
taken the prior probability distributions for /star and tit to be fiat in the logarithm, i.e. 
dPprior/dx oc 1/x. The best fit luminosity function, and constraints on /star and tit/tn are 
shown in the right hand panels of Figure 2 of this supplement for the fiducial model (repeating 
results presented in Figure 1 of the original Letter to Nature) . Marginalizing over /star in this 
distribution provides a probability distribution for tn/tu which is then used in the evaluation 
of limits on the minimum virial temperature Tyir of high-redshift galaxies. 
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3. Scatter in the effective optical depth along lines of sight toweirds different 

quasEirs at 2; > 6 

At a fixed redshift the local value of Lya optical depth (ti) scales as 

T, oc p^iJ-\ (2) 

where pnr is the density of neutral hydrogen and J is the ionizing background intensity. 
However the effective optical depth {tcs) is measured across some path-length through the 
intergalactic medium (IGM) and must be computed from t^s — — logF where F = (e""^') 
is the transmission averaged over the path-length. While the ionizing background will be 
sensitive to cosmic variance in density contrast calculated in 3-dimensional volumes (A) on 
the scale of the mean-free-path for ionizing photons i?mfp, the absorption spectra of quasars 
probe only a narrow cyhnder (skewer) through these volumes. Note that although ti varies 
with z due to the expansion of the universe and the growth factor of density perturbations, 
we focus here on the scatter in t^s at a fixed z. 

We calculate the predicted scatter in r^s using a combination of cosmic variance among 
different regions and Poisson noise in the number of ionizing sources. Our calculation is 
based on an extension of a detailed model that we pubhshed in Wyithe & Loeb (2006). A 
full description of the methodology of the model is given in that paper. Here we give a 
summary, and point out the various additions and improvements. Fluctuations in both the 
ionizing background intensity J about the mean value J as well as inhomogeneities in the 
gas density, contribute to fluctuations in Tefj. The ionizing background is smoothed on the 
scale -Rmfp, since any point in the IGM sees sources within a volume Vmfp = (47r/3)i?^fp. The 
model combines four major sources of fluctuations in Teff : 



(i) Poisson noise in the number of sources contained within a volume V^aifp 

The scatter in the radiation field due to the finite numbers of galaxies is computed by 
weighting the variance in the numl:)er of sources hy luminosity, i.e. 



_V ^2 I_ 

Xw^BCn^a. dMl350L{Ml350)Vmfy^{Ml350) 

where L is the luminosity, and Ali35o,max is the absolute AB magnitude at rest-frame 
1350Acorresponding to the minimum virial temperature Tmm- We define fluctuations in 
J due to Poisson noise to be Sj^p = {J — J)/ J; 

(ii) Fluctuations in the value of J within volumes of Vmip resulting from delayed or enhanced 
structure formation in underdense or overdense regions 
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Large-scale inhomogeneity in the cosmic density field leads to structure-formation that 
is enhanced in over-dense regions and delayed in under-dense regions. The resulting cosmic 
variance in the redshift at which a critical value of collapsed fraction (Fcoi) is reached within 
regions of size i?mfp may be calculated as 

{Szl)^ = ^^{l + z), (4) 

where cr(i?mfp) is the r.m.s. amphtude of the linear density field smoothed over spheres of 
radius i?mfp, and Sc is a critical overdensity for collapse [oc (l + z) at high redshift]. Since the 
ionizing background within a given region depends on the value of the collapsed fraction, we 
can relate {Sz^^)^ to the cosmic variance in redshift where a critical value of J is obtained. 
Wc would like to convert this scatter in redshift to a scatter in the optical-depth to Lya 
absorption at a given redshift. In a smooth (i.e. non-clumpy) IGM the fractional change in 
the value of the optical depth relative to the average, given a density contrast (A^), and a 
delay {Sz) in the redshift where the critical collapse fraction is reached is 



where Tav is the average value of optical depth at z. Within our fiducial model the production 
rate of ionizing photons within a volume V^fp, and hence J, is taken to be proportional 
to the time-derivative of the fraction of baryons which collapsed into galaxy halos (Fcoi). 
There are two contributions to the derivative dJ/dz: (a) the evolution of the star formation 
rate; and (b) the evolution of the co-moving mean-free path. Fan et al. (2006) present 
measurements of -Rmfp at redshifts between 4.5 and 5.65, from which we find R^fp oc (l + z)"' 
with 7 = —1.75 ± 0.75. Hence in our fiducial model we have 



d\nJ _ f d\n{dFcoi/dt) 7 



dz \ dz {1 + z) 

and fluctuations in J due to cosmic variance of 

J— J /dlii{dFcoi/dt) 7 



Sz, (6) 



JJ,cv 



J V dz (1 + ^) 



Sz. (7) 



{in) Fluctuations in Lya transmission due to fluctuations in the density contrast smoothed 
on the scale i?mfp- 

To estimate the typical fluctuations in the effective optical depth in a smooth IGM we 
must compute the variance in overdensity among lines of sight through the density field of 
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length i^mfp in addition to fluctuations in J. We calculate the power-spectrum of fluctuations 
in cylinders of length L — 2Rjnfp and radius R 



where Ph{k) as a function of wavenumber k is the linear baryonic power-spectrum, which 
may be approximated as Pb(^) — + k'^Rf)''^ in terms of the cold dark matter power- 

spectrum P{k) and the filtering scale Rf for the associated reionization history The variance 
in Az on a scale L with wavenumber /cl follows from 



The radius of the cylinder R is set by the size of the quasar emission region, which is 
much smaller than the filtering scale. The fiuctuations in the line-of-sight density field are 
= {A^ — 1) = D5z, where D is the growth factor at z and (with variance (Jl,z) is the 
overdensity smoothed in thin cylinders of length 2R^fp; 

{iv) Fluctuations in the transmission due to small-scale density fluctuations along the line- 
of-sight through a region with constant J 

In addition to the smoothed density field averaged on a large scale, fluctuations in the 
transmission arc sensitive to small scale structure along the line of sight. Inclusion of these 
non-linear effects requires a numerical simulation. The fluctuations in effective optical depth 
ta about a mean 7^ have been computed by Lidz et al. (2006) based on the transmission 
power-spectrum, giving ta = 2.5lo;32) 2.61q;33, ^-Sloisg ^nd 4.0to;48 respectively, in regions of 
co-moving length L — 55Mpc. We extrapolate from these results to compute the fluctuations 
in transmission through thin cylinders of length 2i?infp. The transmission [F) is related to 
the effective optical depth through F — e~'^''^ . Hence, the lower and upper fluctuations in 
transmission (a^) are related to the upper and lower fluctuations in optical depth = T^6r) 
through the relations ap = e~'^'=^{l — e~'^^) and ap = e~'^'=^{e'^^ — 1). On scales larger than 
~ lOMpc the transmission power-spectrum is close to white noise, and ap/F oc L~^/^. To 
estimate fluctuations in ta within thin cylindrical regions of length 2i?infp, we therefore 
multiply the simulated transmission fluctuations (ap/F) by \/L/ (2i?mfp) (treating upper 
and lower limits separately). 

We assume that each line-of-sight samples a probability distribution P(A) of density 
contrasts A relative to the mean density (as described by Miralda-Escude et al. 2000, ApJ, 
530, 1) on scales <^ -Rmfp- Deflning the local optical depth to be t; = AA^, one can flnd an 
expression ta = 'rA(^) for ta given A and P{A) (Songaila & Cowie, 2002, Astron. J., 123, 
2183; Lidz et al. 2006). On the scale -Rmfp, we modulate A based on the fluctuations in density 






(9) 
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and radiation, A' — A{1 + S^f / (1 + 5j;p + 5j,cv), and hence find the modified = '^a(^')- ^ 
distribution of t^s that includes fluctuations in J is then found via Monte-Carlo realizations 
of the values for ta, (^j,p, (^j,cv and ^a- 

Points in the IGM which are separated by 2i?infp have ionizing backgrounds that are 
dominated by different sets of sources. Since the two-point correlation function is much 
smaller than unity on the scale 2_Rjnfp (tens of Mpc), we consider regions of this scale to 
be independent and assume that the transmission power-spectrum remains close to white 
noise in the presence of a radiation field that fluctuates on scales larger than 2i?infp. The 
number of independent regions in the ionizing background within a rcdshift interval 5z is 
A^mfp = 5z{cdt/dz)/{2R^{p), where c is the speed of light. We flnd A^mfp ~ (1.45+o;4o), 
(1.35l2:I^), (1.9lJ:^) and (2.7l?:|]) at z = 5.25, 5.45, 5.65 and 5.85 respectively. To relate 
modeled fluctuations in regions of radius i?mfp to observations made in redshift bins, we 
divide the modeled fluctuations in transmission through regions of volume T4ifp by ^/N^. 
Finally, to facilitate comparison of modeled distributions with separate upper and lower 
fluctuations (corresponding to an asymmetric distribution) and the quoted variance in the 
observed values, we take the root- mean-square of the upper and lower fluctuations in Sr- 

Note that our constraints assume that galaxies dominate J. On the other hand, quasars 
arc more luminous and rare in comparison with galaxies and may therefore contribute to the 
fluctuations in the ionizing background. However, the known population of quasars cannot 
contribute to the scatter among regions of size Vmfp, since there is only ~ 1 quasar as bright 
as the SDSS quasars [or the faintest known X-ray quasar (Bargcr ct al. 2003, ApJ, 584, 
L61)] at 2; ~ 6 per 10^ (or 100) regions of radius -Rmfp ~ 20Mpc. Moreover, the contribution 
from quasars to reionization is thought to be much smaller than is required to reionize the 
universe at 2; ~ 5.5-6. Even if there was one dominant quasar per Knfp, it would need to 
contribute > 10% of the ionizing background in order to influence our results. 



4. A-posteriori constraints on Tmin 

The model described above predicts the fluctuations in the large scale optical depth 
given a set of parameters including T^m, tit, 7, and -Rmfp- By comparing this model to the 
data we flnd constraints on Tmin using prior information on the other parameters. To compute 
the distribution dP/dT^^^, we flrst flnd the likelihood for Tmin, which may be computed as 



T ■ oc 

min 



/■ dP r , dP f , dP ^ ^ , 

J o-Kmfp^^-y J d-f— J dtit—Ls^{Rmfpn,tit,T^in)- (10) 
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Here the a-posteriori distribution for lifetime is computed from ^ oc J c^/star^^^-^^, and 
-^^Y is determined from Fan et al. (2006). The prior probabihty distributions for /star and 
t\i are taken to be flat in the logarithm. The likelihood Lj^(i?iiifp. 7, tit, Tmm) was computed 
from the observed value of the variance in 5^. The imcertaintics on 5r were computed 
assuming Gaussian scatter in t^q and a sample size of 19 at each redshift (corresponding to 
19 lines-of-sight). Wc find the a-posteriori distribution for T^\^ using a flat prior probability 
distribution in log(Tinin), i.e. ji^, oc L 



5. Constraints on Tmm including merger driven stEir formation 

For simplicity our fiducial analysis assumed that the star formation rate is approxi- 
mately proportional to the derivative of the collapsed fraction of baryons in halos with a 
virial temperature Ty^j. > T-^^in- However, starburst episodes are believed to be triggered by 
galaxy mergers. We have therefore repeated our calculation of the constraint on T^m with a 
prescription for star formation that is based on galaxy mergers rather than the derivative of 
the collapsed fraction. The star formation model combines the contribution from mergers of 
halos with an initial Tyir < ^min (that go above Tmin after the merger) with the contribution 
from mergers of galaxies with an initial Tyu > Tmm. Thus, gas that has not previously cooled 
and undergone an episode of star formation in a galaxy forms stars with efficiency /star as it 
cools following a merger to form a galaxy more massive than Mmin = M{T^i^). In addition, 
gas in a galaxy of mass M > Mmin undergoes an additional starburst if the galaxy merges 
with a second galaxy whose mass is larger than M/2. The model assumes one previous burst 
of star formation in the gas involved in the major merger so that a fraction 1 — /star (where 
/star is the star formation efficiency) of the gas is available for the new burst. We assume 
/star = 0.1, though the results are quite insensitive to this parameter. The resulting star 
formation rate may be written 



min(M,Mniin) at 

, » , ^ u i V mrg 



jp O poo /'mm(M,lvi^i^) 

^{z) = ^ / dM dAM 

dt iim iM^i„/2 imax{0,M^in-M) dAMdt 



^^fsi.r (e(Mmm - M)M + AM) 



+^ /- m r dAM 



dt dn 
^ dzdM 



^^m Jm^^ Jm/2 dAMdt 

X/sta.(l - /sta.) [e(M - Mmi„)M + e(AM - Mmin) AM] , (H) 

where Qra and flh are the fractions of the critical cosmological density comprised of mass and 
baryons respectively. Here -^{z) is the Press- Schechter mass function of dark-matter halos 



at redshift z, and \m the number of mergers per unit time of halos of mass AM with 

halos of mass M (forming new halos of mass Mi = M -|- AM) at redshift z. The function 
Q{x) is the Heaviside step function. The first term accounts for star formation when gas 
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Fig. 2. — Constraints on the star-formation efRciency (/star), duty-cycle (tit/tH) and minimum virial 
temperature (Tmin) of galaxies at z 5.5 — 6. The constraints were placed assuming the merger driven 
model for star formation. Panel a: The luminosity function inferred from the UDF, together with the 
best-fit solution for a luminosity function model based on the Sheth-Tormen mass function (with luminosity 
proportional to halo mass). The upper axis shows the halo mass corresponding to A^i35o for the best fit. 
Panel b: Contours (at 64, 26, and 14% of the maximmn likelihood) of the joint a-posteriori probability 
distribution [(PP/d{\og f star) d{\og tit)]. For comparison, the grey contours show the virial temperature of 
halos [labeled using log]^Q(Tvir)] corresponding the lowest luminosity bin in the UDF. Panel c: The joint 
probability distributions for log (tit/in) and log (Tmin), which include the constraint on tit/<H from the galaxy 
luminosity function. The contours represent values at 64, 26 and 14% of the maximum likelihood, and are 
shown for observations at z ~ 5.45 (blue contours) and z ~ 5.65 (red contours). This distribution includes 
constraints on t\t from the galaxy luminosity function. However assuming the most likely value of i?mfp, we 
find that T„iin > lO^K at z < 5.5 in all cases where the duty-cycle of the UV-bright phase of galaxies is greater 
than 1% (corresponding to the lifetime of massive stars). Panel d: Differential probability distributions 
for log (Tniin) at z ~ 5.25 (dashed blue lines), z ^ 5.45 (blue lines), z ^ 5.65 (red lines) and z ~ 5.85 
(dashed red lines). We also show combined constraints using the lowest three redshift bins (black lines). For 
comparison, the upper axis shows the corresponding masses at z = 5.5. In this case the solid and dotted 
lines correspond to differential and cumulative distributions. Throughout this work, we adopt the latest 
values for the cosmological parameters as inferred from the Wilkinson Microwave Anisotropy Probe data. 
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first assembles into a halo with M > Mmin- The second term corresponds to the secondary 
starbursts (in gas that has already undergone a burst of star formation) that follow major 
mergers. 

This alternative estimate of the star formation rate leads to results that are shown in 
the left panels of Figure 2 of this supplement. The preferred value is Tmin ~ 10^-^''=°-^K, with 
Tlnin > lO^K at the > 99% confidence level. 



6. Constraints on T^am including supernovae feedback 

In addition to the suppression of star formation in dwarf galaxies due to a reionized 
Universe, there may also be suppression of star formation within low mass galaxies as a 
result of supernovae driven-winds which expel gas from their shallow potential wells. In the 

local universe, this suppression is apparent at observed circular velocities 14bs 178 km s~^. 
A recent study has shown that in a large sample of local galaxies, the ratio e = M^/Mhaio 
(where and Mhaio are the total mass of stars and the dark matter halo respectively) scales 
as 

e oc {M^,,jMUf' for M < 

oc 1 otherwise (12) 

The star formation efficiency is proportional to e. 

To explore the effect of supernovae feedback we have repeated our calculation of the 
luminosity function with the following modification. We allow the critical mass M^^^^ at 
z ~ 5.5 to be a free parameter, and fit the luminosity function in equation (1) of this 
supplement, with the derivative ^^^^^^^ modified from the fiducial case because L oc eMhaio- 
Since feedback is thought to operate over a fixed fraction of a dynamical time, we assume the 
starburst lifetime to be constant among galaxies of different mass. To avoid discontinuities 
in the derivative ^^^^^ , we smooth e(M) using a log-normal window function of width 
0.25dex. The best-fit luminosity function model that includes supernovae feedback is shown 
in Figure 3 of this supplement. We find that allowance for supernovae feedback suggests high 
redshift galaxies to be located in more massive halos, with long duty-cycles tit/tu > 0.5. 

Estimates of the scatter in the ionizing background due to Poisson noise in the number of 
galaxies and due to fluctuations in the star formation rate resulting from delayed or enhanced 
structure formation must also be modifled to account for supernovae feedback if it operates 
at high redshift. To explore the effect of supernovae feedback on fluctuations in the ionizing 
background and hence on conclusions regarding the suppression of dwarf galaxy formation 
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Fig. 3. — Constraints on the star-formation efSciency (/star), duty-cycle {tit/tn) and minimum virial 
temperature (T^i^) of galaxies at z ~ 5.5 — 6. The constraints were placed assuming the fiducial model for 
star formation, with suppression in dwarf galaxies due to supernovae feedback. The lines and symbols have 
the same meaning as in Figure 2. 



in the reionized IGM, we have repeated the comparison of our modeled fluctuations with 
the observed fluctuations in Tes, including the following modifications. First, in calculating 
the Poisson noise we include a galaxy luminosity function that takes account of supernovae 
feedback in both the numerator and denominator of equation (5). Second, our fiducial 
model assumes that the derivative of the dark-matter collapsed fraction multiplied by the 
star formation efficiency (/star) and baryon fraction, approximates the star formation rate. 
We replace the derivative of collapsed fraction in this formalism with 

—P- = — / dMeM—, 13 
dt dt Jm . dM' ^ ^ 

-"J mm 

where e is computed using equation (12) of this supplement. The results are shown in 
Figure 3 of this supplement. The most likely value is Tmin ~ lO^'^-o ^K in agreement with 
our previous calculation, while Tmin ^ IC^K is ruled out at the > 95% level. 



